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Abstract 



I We investigate a connection between a property of the distribution and a conserved 

I quantity for the reversible cellular automaton derived from a discrete-time quantum 

■ walk in one dimension. As a corollary, we give a detailed information of the quantum 

_ , walk. 
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1 Introduction 



Recently considerable work lias been done on quantum walks (QWs), wliich 
may be useful in constructing efficient quantum algoritlims. Kempe [1], Tre- 
genna et al. [2], Ambainis [3] liave reviewed the field. For more general setting 
in connection with quantum cellular automata, see Meyer [4]. Two classes of 
QWs have been studied: one is a discrete-time case, the other is a continuous- 
time case. Here we focus on discrete-time QWs in one dimension. The QW 
can be considered as a quantum analog of the classical random walk. However 
there are some differences between them. For example, a QW starting from 
the origin at time n = moves away faster than the classical random walk. 
The standard deviation a{n) of the probability distribution for the QW in- 
creases with time according to o"(n) ~ n, while for the classical random walk, 
cr(n) ~ ^/n. 
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In this Letter we consider some properties on a one-dimensional reversible 
cellular automaton (RCA) derived from a QW on a line. We present necessary 
and sufficient conditions on the initial state for some conserved quantities of 
the RCA. One is the expectation of the distribution (the 1st moment), and the 
other is the squared norm of the distribution (the 0th moment). The former 
corresponds to the symmetry of distribution, because our result (see Theorem 
2) implies that the distribTition is symmetry for any time step if and only 
if its expectation becomes zero for any time step. We should note that our 
RCA is different from a quantum cellular automaton studied by Grossing and 
Zeilinger [5,6]. In [6], they found a conservation law connects the strength 
of the mixing of locally interacting states and the periodicity of the global 
structures. Their conservation law also differs from our conserved quantities. 
Concerning a recent review on quantum cellular automata, see Aoun and Tarifi 
[7], and Schumacher and Wernerfor [8], for examples. It is well known that 
cellular automata are models based on simple rules which upon deterministic 
time evolution exhibit various complex behavior [9]. So our investigation on 
the RCA might be applicable to the analysis on the behavior. Moreover, by 
applying our results to the QW case, we obtain a more detailed information 
on QWs. 

The rest of the Letter is organized as follows. Section 2 gives the definitions 
of the QW and the RCA. In Section 3, we study a necessary and sufficient 
condition on the initial state of the symmetry of distribution (in general, non- 
probability distribution) of the RCA. Section 4 treats a necessary and sufficient 
condition on the initial state for a conserved quantity of the RCA. In the final 
section, we apply our results to the QW. 



2 Definitions of the QW and the RCA 



In this section first we present a definition of the QW. The time evolution of 
the QW is given by the following unitary matrix: 



U 



a h 
c d 



where a,b,c,d E C and C is the set of complex numbers. The unitarity of 
U gives |ap + |cp = \b\^ + \d\^ ^ I, ac + bd = 0, c ^ -Ab and d = Aa, 
where ^ is a complex conjugate of 2; G C and A = det U = ad — be. The QW 
is a quantum generalization of the classical random walk with an additional 
degree of freedom called the chirality. The chirality takes values left and right, 
and means the direction of the motion of the particle. At each time step, if 
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the particle has the left chirahty, it moves one step to the left, and if it has 
the right chirality, it moves one step to the right. 



The evolution is defined by the following way. The unitary matrix U acts on 
two chirality states \L) and \R): \L) a\L) + c\R) and \R) b\L) +d\R), 
where L and R refer to the right and left chirality state respectively. In fact, 
put \L) = *[1,0] and \R) = *[0, 1] where t means the transposed operator. 
So we have U\L) = a\L) + c\R) and U\R) = b\L) + d\R). More in detail, at 
any time n, a 2-vector e at the location A; e Z, where Z is the set of 
integers, gives the amplitude of the probability of finding the particle at k: the 
probability is given by the square of the absolute value of the vector at A;. If 
l^jt(n)) define the amplitude at time n at the location k where 



n 



with the chirality being left (upper component) or right (lower component), 
then the dynamics for \^}.{n)) in the QW is described by the following trans- 
formation: 



\^k{n + 1)) = \L){L\U\^k+i{n)) + \R){R\U\^k-i{n)) (1) 

for any time n G Z+, where Z+ is the set of non-negative integers. It is known 
that the evolution of the QW is highly sensitive to the initial qubit state (see 
Refs. [10,11,12,13], for examples). Therefore we define the set of initial qubit 
states as follows: 



$ = {(^ = * [^0^(0), *^(0)] e : 1*0^(0)1^ + K(0)r = 1} 



Here we introduce some typical QWs. The simplest and well-studied QW is 
the Hadamard walk whose unitary matrix U is defined by 



1 1 
1 -1 

The dynamics of this walk corresponds to that of the symmetric random walk 
in the classical case. However, it is noted that the symmetry of the distribu- 
tion for the Hadamard walk depends heavily on the initial qubit state, see 
Konno et al. [10] (the Hadamard case), Konno [11,12,13] (general case). An- 
other example is an extension of the Hadamard walk in the following way (see 
Romanelli et al. [14], for instance): 



H 



1 

71 



3 



cos 6 sin 6 
sin 9 — cos 9 



where 9 G [0, 27r). In this Letter we consider mainly the case of ^ G [0, 7r/2) for 
the sake of simphcity. Remark that i7(7r/4) = H. Now Eq. (1) can be written 



as 



Furthermore, we can apply the above equations to each other to uncouple the 
chirahty components. Resultantly, both ^i^(n) and ^i^(n) satisfy the following 
same partial difference equation: 



Xk{n + 2) = aXk+i{n + 1) + dXk-i{n + 1) - {ad - bc)Xk{n) (2) 

The above argument appeared in Knight et al. [15,16] (see also page 279 in 
Gudder [17]). In other words, Eq. (2) imphes a dynamical independence of 
the evolution of the two chiralities L and R. Thus there are two essentially 
independent walks, coupled only by the first two steps. After that the two 
walks can behave independently each other. In this Letter, we call the cellular 
automaton defined by Eq. (2) reversible cellular automaton (RCA) , since Eq. 
(2) imphes that Xk{ri) can be also determined by Xk{n+ 1) and Xk{n + 2). 
In the terminology of cellular automata, the state space at each lattice site is 
usually taken to be finite, however, continuous case can also be chosen (see 
pp. 155-160 in [9], for instance). A special case of a = d = 1 and b = c = 
is equivalent to a discretized two-dimensional wave equation derived from 
numerical modeling of cosmic strings, see Ref. [18]. 

Here we give an important relation between the QCA and the QW. If we put 
(Xo(0),X_i(l),Xi(l)) = (^^(0),a^^(0) +6^^(0),0) on the initial state of 
Eq. (2), then Xk{n) — ^^{n), that is, Xk{n) represents the left chirality of 
the QW. Similarly, if (Xo(0), X_i(l), Xi(l)) = (*^(0), 0, c^^(O) + (i^^(O)) is 
chosen, then Xk{n) — ^-^(n), i.e., Xk{n) represents the right chirality of the 
QW. 

In this Letter we consider a more general setting concerning the initial state. 
So we define the set of initial states for the RCA: 

$ = (Xo(0),X_i(l),Xi(l)) = (a, A 7) : a,A7 e C} 

To distinguish our initial state for the RCA considered here from that for the 
QW in our previous papers [10,11,12,13], we put "~" in $ etc. Prom now on 



4 



we consider only the RCA defined by H{9): 



Xk{n + 2) = cos9[Xk+i{n + 1) - Xk-i{n + 1)] + Xk{n) (3) 

Here we define a distribution (in general, non-probability distribution) of the 
RCA X{n) at time n e Z+ by {\Xk{n)\'^ :keZ}. 

McGuigan [19] investigated some classes of quantum cellular automata. In his 
setting, Eq. (3) belongs to a class of fermionic quantum cellular automata 
whose update equation is given by 



X,{n + 2) = /(Xfc+i(n + 1), Xfc(n + l),Xk-i{n + 1)) + Xfc(n) 

Our case is f{x,y,z) = (x — z) cos9. By using Eq. (3), Romanelli et al. 
[14] analyzed in detail discrctc-timc one-dimensional QWs by separating the 
quantum evolution into Markovian and interference terms. 



3 Symmetry of distribution of the RCA 

First we present the following useful lemma to prove Theorem 2. 
Lemma 1 (i) Suppose that initial state is 

^= (Xo(0),X_i(l),Xi(l)) = {a,f3,-f3) 

where a, (3 e C. Then, for any k,n E Z^, 

Xkin) = (-l)"X_fe(n) (4) 

(ii) Suppose that initial state is 

(p^{Xo{0),X_,{l),X,{l)) = {0,P,e'^P) 

where /3, ^ e R and R is the set of real numbers. Then, for any k,n e 

Xfc(n) = (-ir+v^x:;H (5) 

The above lemma directly comes from the symmetry: if Xk{n) is a solution 
of the RCA, then so is xl^\n) = (— l)"X_fc(n), and also so is, for any ^, 
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X^\n) — (— l)"+^e'^X_fc(n). To state Theorem 2, we introduce the following 
three classes: 

$x = {(^e$:/3 + 7 = (6) 

or = |7|(> 0), and q; = (7) 
or \P\ = |7|(> 0), a 7^ 0, and 0p + 0^- 29^ = n (mod 2n) } (8) 
= e $ : \Xk{n)\ = for any n E Z+ and /c e Z} 

{oo 
<^ G $ : ^ A;|Xfc(n)|^ = for anyn e Z^ 
fe=— oo 

where 9z is the argument of 2; G C with 2; 7^ 0. It is noted that if G $s, then 
the distribution of X(n) is symmetric for any time step n G Z_|_. 

Theorem 2 For ^/le RCA defined by Eq. (3), we have 

Proof. First we see that the definitions of $s and $0 give immediately 

C $0 (9) 

Next we prove $_l C $5. The proof is divided into the following three cases. 
Case 1. If an initial state satisfies Eq. (6), then we can obtain the following 
relation from Lemma 1 (i): 

\Xk{n)\^\X_k{n)\ (10) 

Case 2. Assume that an initial state satisfies Eq. (7). Then Lemma 1 (ii) 
implies Eq. (10). Case 3. Assume that an initial state satisfies Eq. (8). Let (pi ~ 
(|a|e*^", |/5|e^^^ \j3\e'^'^-^^+'^^'^^). So we have ^2 = e-^^«(^i = (|a|, \j3\e'^^^-^"\ 
|^|gi{^-e^+e„))_ We put e = Op-Oa,, then (^2 = (|a|. |/^|e*(^-^)). Therefore, 

real and imaginary parts of (^2 become 'Rc{(p2) = (|a|, cos^^, — cos ^) 
and Im((^2) = (0, \I3\ smO, \ j3\ sin^^), respectively. Here if we take Re((^2) as an 
initial state of the RCA, then we can see that distribution of X(n) is symmetric 
by using Lemma 1 (i). Similarly, using Lemma 1 (ii), the initial state Im(<^2) 
also gives symmetric distribution. In order to confirm (p2 = Re((^2)+^Ini((^2) £ 
$5, we can use the linearity of Eq. (3) and the fact that if the initial state 
{a,f3,'j) is real, then the solution {Xk{n)} is also real. From (^2 G ^s, we can 
get (pi G $s. Therefore for any ip G we have |X;i.(n)| = |X_fe(n)| for any 
A; G Z and n G Z+. Then we conclude 
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Finally, a direct computation gives 



(11) 



m(i) = i7r - i/?r, m(2) = 2 cos^ ^(i7r - \m 

m(3) = ^(3cos^2^ + 2cos2^ + l)(|7|2- l/^n 

-^sin6' sin 26' {q;(^+7) + + 7)} 

00 

where ■m{n)= ^ A;|Xfe(n)|^. The above equations imply that if <^ e $0 then 

A;=— cx) 

e $_L. So we have 



^0 C $x (12) 

Combining Eqs. (9), (11) with Eq. (12) implies I>_l C $s C $0 C so the 
proof of Theorem 2 is complete. 



4 Conserved quantity for the RCA 

From now on we assume that < ^ < 7r/2 for simplicity. First we let 



k=—oo 

Here we consider a necessary and sufficient condition on the initial state for a 
conservation law: ||X(n)|| — c{n >0) with c > 0. To do so, we introduce the 
following two subsets of 



$*(c) = {(^ e $ : \a\^ = c, (13) 

m'+b\'^c, (14) 

/37_+/37 = 0, (15) 

q;(/3 - 7) + a(/3 - 7) = 2c cos e } (16) 

and 

$(c) = 1^ e $ : ||X(n)||^ = c for any n e Z+| 
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where $(c) stands for the set of initial states statisfying that ||X(n)|| becomes 
a conserved quantity. When c = 0, it is easily shown that $(0) = $*(0) = 
|(^ e $ : q; = /? = 7 = o| . So we assume that c > 0. Then we have 

Theorem 3 For any c> 0, 

Mc) = He) 



The above theorem implies that $*(c) gives a necessary and sufficient condition 
we want to know. 

Proof. Case 1: $(c) C $*(c). A little algebra reveals 



\\xm\' = H', (17) 
ii^(i)ir = i/3r+i7r, _ m 

||X(2)||2 = + 2cos'^(|/5|2 + - cos" 6 {(3^ + [3 

-cose{a(P-^) + a{P-'y)}, (19) 
||X(3)|p = 2cos2^|a|2 + {2cos^^ + (1 -2cos2^)2} + \^\^) 
+2cos2^(l - 2cos2^)(/37 + ^7) 

+ cos^(l -3cos^^){a(;5-7) + a(/3-7)} (20) 
Prom Eqs. (17) - (19) and < ^ < 7r/2, we have 



a{P - 7) + a(/3 - 7) = cos e{2c - {(5^ + (3^)} (21) 

Combining Eq. (21) with Eqs. (17), (18), (20) implies sin6' cos 6'(/?7 + ^7) = 
0, so we have I3q + [3^ = 0, since < ^ < 7r/2. Then Eq. (21) becomes 
a{P — 7) + a{P — 7) = 2c cos 9. So we have the desired result. 

Case 2: $*(c) C $(c). Define a spacial discrete Fourier transformation of Xfc(n) 

by 

fc=— 00 

where ^ G R. Prom Eq. (3), we have 



where Xo{^) = a and Xi{^) = e~*^/9 + 6*^7. We define ip = ip{^) G R by 

A+(^) = e'^'^, that is, cosy? = y^l — cos^ 6 sin^ ^ and sin 93 = — cos 6* sin ^. So 
A_(0 = -e-''^. Then we obatin 
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where A±(^) = —i cos 6*8111^ ± ^JY^-cos^'o'^j^ and 
Define 





In particular, we put ||/||* = (/I/)- Then 

=-oo 



fe=— oo 

imphes | \X{n) \ \ = \ |-^n| |*- In this setting, we have the following result: for any 
n e Z+, 

||X(n)||^ = \\A\\l + + {-iy'{{e''"^A\e~''"^B) + {e-'""^ B\e'''^ A)} 
A direct computation gives 

u\\i = m\i 

/ \ 2 

/ n -J, . I 1 — sm6\ 
-(/57 + /37) 



cos^ 



Moreover we have 



I |2 /• cos(2(n - l)ip) 
= \a\ / 5 

COS^ (/? 







27r 27r 

+ (,j^_,^)J ^ eosesin((2n-l)^) . sinesin((2n - 1)^) 

J cos^ CP 7 cos^ (fi 



COS^ </? J cos^ ip 
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+ a 



27r 27r 

I . cosesm((2n-l)c,) . sinesin((2n - 1)^) 

1^ 7 COS"' </? J COS"' 

J COS^ iO 
^ 

I J cos{2i)cos{2nip) _ . J sm{2^) cos{2n(p) 

I y COS^ (f J COS^ </? 

{27r 27r 
/• cos(20 cos(2ny) ^ • / sm(20 cos(2ny) 
J COS^ </? J COS^ (f 



COS^ </? J COS^ (p 

27r 

cos(2n(/9) 

COS^ iO 
^ 



for any n e Z+. Using a change of variable in the integration and noting 
r cos^sm{{2n-l)(p) r sm{2^) cos(2ny?) ^ q 

J COS^ <Z) J COS^ 



we have 



1 



2 sine 



-H/3-7)+a(/3-7)} 



1 — sin 9 ' 
cos^ 



/ ^ T> X 1 ~ sin ^ ' 
-(/37 + /37) 



cos 6* 



(-ly 



TT 



e-Tv/2 



cos(2(n — l)a;) 
cos xV cos^ x — sin^ 6 



dx 





, ,- , ,^ M 1 / sinxsin((2n - Ijx 

-{«(/3-7)+a(/3-7)}— - / '\ 

cos (7 J ms.TVCos^ .T — sm 



-(/37 + ^7)< 



e-Tr/2 

cos(2n,T' 



6»— 7r/2 


cosxV 


f 




/ - 


2 


cos^ 






6>-7r/2 









-/ 

6>-7r/2 



COS cos^ x — sin^ Q 
dx 



dx 



dx 



cosx 



cos( 2ii.i 



cos cos^ X — sin^ ^ 



(22) 
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for any n G Z+. The condition "jap = + = c, q;(/3 — 7) + a{P — 7) 
2c cos^, /?7 + ^7 = 0" in $*(c) gives 



,2 , ? K{x) 



\\X{nW = c + ^-^ / y ; ' . , dx 

J , COS XV cos^ X — sm 

0--K/2 

where hn{x) = cos(2(n — l)x) — 2 sinx sin((2n — l)x) — cos{2nx) . On the other 
hand, a httlc algebra reveals that hn{x) = for any n G Z^. Therefore we 
conclude that ||X(n)||^ = c (n = 0, 1, 2, . . .), that is, $*(c) C $(c). The proof 
of Theorem 3 is complete. 

By using Eq. (22) and the Riemann-Lebesgue lemma (see page 462 in Durrett 
[20], for example), we have 



lim \\X(n)"'^ 



2 sin 6* 



1— sin^A ,^ — , fl — sm6\^ 



-MP - 7) + - 7)} w- - (/37 + /37) 



cos^ / V cos^ 



(23) 



The above equation will be used in next section. As a corollary to Theorem 
3, we can show that any RCA with symmetric distribution does not have our 
conserved quantity: 

Corollary 4 For any c> 0, 

$5 n $(c) = 



Proof. We consider the following three cases. Case 1. When /3 + 7 = 0, from 

l/^P + I7P — c, we have P = ^J'c/2e^^i^ and 7 = —\Jc/2 e*^". However the last 
equations and j3'^ + (3^ = contradict each other, since + f3'~f = —c {< 0). 
Case 2. = |7|(> 0), and a = 0" and "lap = c(> 0)" contradict each other. 
Case 3. We assume that = |7|(> 0), a 7^ 0, and^^+^^-2^c« = 7r(mod27r)". 
Prom now on we omit "mod 27r" . Then q;(/3 — 7) + a{P — 7) = 2c cos 9 can be 
rewritten as 



\a\{\P\ cos{9a - 0(3) - I7I cos(^a - 9^)} = c cos^ (24) 
On the other hand, \a\ — y/c, \P\ — \ j\ — ^Jc/2. So Eq. (24) becomes 

COs(^a - Ofi) - COsiOa - 9^) = V2 COS 9 (25) 
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Combining Eq. (25) with 9^ + 9^ — 29a — tt implies 



sin {{9f3-9^)/2) = -cose/V2 



(26) 



We see that P^+P'y = gives 9p—9^ = ±7r/2. So Eq. (26) becomes cos6' = ±1. 
Then we have a contradiction, since we assumed that 0<^<7r/2. So the 
proof of Corollary 4 is complete. 



5 The QW case 



In the last section we return to the QW given by H{9). So we put initial qubit 
state (p = ^[ai,ar] with \ai\'^ + = 1 at the origin for the QW. Here we 
apply our result (Theorem 3) on conservation property of the RCA to the QW. 
First we consider the left chirahty case, that is, a — ai, (3 — cos 9ai-\- sin 9 Ur 
and 7 = 0. Then Eq. (15) holds, since 7 = 0. Eqs. (14) and (16) can be 
rewritten as 



cos^ 9\ai\'^ + sin^ ^Ictrl^ + cos^sin^(a;Q;r + aja^) = c 
2cos^|q;/|^ + sin 9 {aic^ + ajar) — 2c cos 9 



(27) 
(28) 



respectively. By using sin ^ 7^ (since < ^ < tt/2) and Eqs. (13) and (28), 
we have 



aiUr + aiar — 
Combining Eq. (27) with Eq. (29) gives 



(29) 



cos^ 9\ai\ 



sin^ 9\aJ'^ = c 



On the other hand, it follows from I a; I 



c. By Eq. (30) and the last equation, we get 



1/2 



(30) 

1 and Eq. (13) that ja^p = 



(31) 



Moreover \ai\ = \ar\ — l/V^ can be derived from Eq. (31). So combining this 
result and Eq. (29), we conclude that the initial state of the QW given by the 
RCA having our conserved quantity with c = 1/2 can be determined as 



ai 



± 



V2 



(32) 
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where ^ G R. Similarly, we consider the right chirality case, that is, a — 
a,., /? = 0, and 7 = sin 6 ai — cos 6 ar- In this case also, a similar computation 
gives the same conclusions, i.e., Eqs. (31) and (32). 

On the other hand if we choose Eq. (32) as an initial qubit state for the QW 
defined by H{6), then the probability distribution becomes symmetric (see 
Refs. [10,11,12]). Furthermore our result gives an additional information on 
the conserved quantity: 

\\^\n)\\' = \\^^{n)\\'=^- 



for any n > 0. In particular, when 6 = n/A (the Hadamard walk), ii (p — 
^[1/V2,i/V2], then ||^^(n)||2 = ||^^(n)|p = 1/2 for any n > 0. However 
Corollary 4 implies that the conserved quantity (in our case c = 1/2) is not 
compatible with the symmetry of the distribution. This gives an interesting 
result on the QW: any symmetric probability distribution of the QW can not 
be written as the sum of the symmetric distribution for each chirality of the 
QW. 

Finally we consider a general case. The unitary evolution of the QW implies 
||*'^(n)|p + ||^'^(n)|p = 1 for any time step n e Z+. Furthermore Eq. (23) 
gives 



lim 11^ (n) 



'1 + cos^6')|a,p + sin^6'|a 



n-*oo 2 sin ^ 

+ sin 9 cos 9{ai'a^ + aiar) 

—{2 cos^|q;;|^ + sin^(Q;/a^ + oTiar)} 

1 



1 — sin I 
cos 9 



lim ||*V)ir= n . . 
„_ooii ^ 2sin^ 



1 + cos^ 6')|Q;r|^ + sin^ 6'|Q;i|^ 
— sin 9 cos 9{aia^ + oqar) 

—{2 cos^^la^p — sin 9{aia^ + alar)} 



1 -sin^^ 



cos 9 

For example, if we choose 9 G (0, 7i/2) and (p = *[1, 0] (asymmetric case), then 
lim ||^^(n)||2 = l-^G f-,lV lim ||v^«(n)|p = ^ G (o,-) 
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